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Abstract 



Several topics of relevance to low transverse momentum i/j and xi,2 (cc) pro- 
duction in polarized proton-proton collisions are discussed. The leading O(o|) 
contributions to the low px Xi production cross-sections via gg, qg, and qq initial 
states are calculated as well as the corresponding spin-spin asymmetries. We find 
that xi production increases relative to direct ip and \2 production, providing 
up to 25% of the observable e + e~ pairs arising from ip decays in pp collisions at 
i/i = 500 GeV . The spin-dependence of xi production, however, is much smaller 
than for either direct ip or \2 production and so will likely be far less useful than 
either process in probing the polarized gluon structure function of the proton. A 
subset of the O(o|) radiative corrections to \2 production involving initial state 
quarks are also performed and compared to leading order gg —> \2 predictions. 



1 Introduction 

Even before the recent approval [[l] of a program of polarized proton-proton collisions 
at collider energies at RHIC |2|, |3], Q , there had been a rapid increase in the literature 
discussing the prospects for probing the spin-dependence of many QCD processes and 
their sensitivity to various polarized parton distributions in such high energy polarized 
hadron collisions. 

Familiar processes such as jet production || ||], direct photon production || [7], and 
weak boson production M, |9| have all been examined in this context in some detail. 
For a few processes, such as polarized Drell-Yan production [fH] and direct photon 



production [11], the full set of radiative corrections have even been computed. In these 
cases, the next-to-leading (NLO) order corrections retain almost all of the features 
of the leading-order (LO) calculations, including their sensitivity to the sea quark 
and gluon polarizations respectively. On the other hand, heavy quark production has 
recently been studied |12|] and it seems likely that NLO corrections (especially at large 
Pt) will change the spin-structure of 6-quark production rather dramatically, pointing 
out the importance of the comprehensive study of NLO corrections to spin- dependent 
processes which might be measured with the proposed RHIC detectors. 



Quarkonium (especially J/ip) production, both at low transverse momentum [13, 



H] , |l"5| , |16ll and at high px [|T7|, |I8|1 , has also been examined for its possible spin- 



dependence and sensitivity to the polarized gluon content of the proton. One intriguing 

suggestion, due to Cortes and Pire |T3|] (hereafter CP), is to consider low px X2i.cc) 

production where the dominant lowest-order subprocess would be gg — > \2- To this 

order, the observable polarization asymmetry, 

_ K++)-a(+-)) 
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(where ± refers to the proton helicity) in such quantities as da/dy(y = 0) (or da/dxp) 
is directly proportional to (AG(x)/G(x)) 2 times a partonic level spin-spin asymmetry 
a 1,1,. Here AG(x) = G + (x) — G-(x) is the longitudinally polarized gluon density, where 
+/— refers to a gluon with its helicity in the same/opposite direction as that of the 
proton. The partonic level asymmetries for X2/X0 production have been calculated in 



the context of potential models |16] and are maximally large, clll = — 1/+ 1 for X2/X0 
respectively, as expected from simple angular momentum conservation arguments. The 
study of the angular distribution of the X2 -^ i>l decay (which also provides the clean 
signal) can simultaneously measure the value of clll O as a further check. Similar 
experiments using xo production are hampered by the very small radiative branching 
ratio to ^7 final states. The spin-dependence of the other dominant if) production 
process at low pt, the so-called 'color-bleaching' mechanism, gg —>■ gi/j, has also been 
examined p6| . 

In the context of the RHIC spin physics program, the PHENIX detector [Qj, with 
its excellent electron and photon detection, will likely be able to resolve ip and Xi,i 
states very effectively so that further study of ip and xi,i production mechanisms at 
low pt in proton-proton collisions (polarized or not) is quite timely. (We note that the 
CDF Collaboration has presented promising preliminary results on such an analysis fl~9 



demonstrating the ability to separate (at some level) direct ip from radiative Xi,2 decay 
in a high energy collider environment at large transverse momentum.) 

Questions can easily arise, however, concerning the reliability of the theoretical 
understanding of x(QQ) production mechanisms in hadronic collisions. Issues such 
as the applicability of potential model analyses of quarkonium production in hadronic 
collisions at low values of Q 2 = M 2 (especially for charmonium) and, as mentioned 



above, on the role of higher order QCD corrections are obvious ones. Both low trans- 
verse momentum inclusive ip and T production |2l], ^TJ and high px ip production 



£| have been well described by such models which do, however, require 
relatively large K factors, namely K « 2, not unlike standard lowest-order jet, direct 
photon and heavy quark analyses. Recent calculations of radiative corrections to low 
Pt 1 Sq(QQ) (i.e., r/g) production in hadron collisions [p6fl , albeit for toponium pro- 



duction, confirm the presence of relatively large K-factors. These corrections would 
presumably be even larger for r\ c states because of the larger value of as at the lower 
charmonium mass scale. The 0(a s ) corrections to the inverse process, i.e. \2 -^ 99 
have been known for some time |27] and were found to be quite large also indicating 
the possible importance of higher-order radiative corrections in P-wave quarkonium 
processes. In addition, data on the relative yields of ip and xi,2 states in pp collisions 



also seem to be compatible with such potential models as well [^8 . 

The presence of large multiplicative i^-factors is, by itself, not a strong argument 
against the use of the CP method as their analysis relies simply on the 2 — > 1 kine- 
matic structure of the gg — > \2 subprocess. Purely virtual corrections to this process 
(retaining the 2 — > 1 kinematics) are unlikely to change the LO predictions for the spin- 
asymmetry which is dictated by helicity conservation. The corrections arising from the 
2 — > 2, gg — > \i9 process at the same order (which are required to cancel the infrared 
divergences found in the virtual diagrams) can, however, change the kinematic struc- 
ture of X2 production and hence the underlying spin-spin asymmetry. Perhaps just as 
importantly, the NLO 0(a s ) processes qg — > q\2 and qq — > g\2, present at the same 
order, introduce a new dependence on the quark and antiquark distributions so that 
even if the contribution of quark/ antiquark initiated processes is relatively small in the 



unpolarized cross-section, it could make a significant impact on the spin-dependence 
depending on the relative sizes of AG(x) and Aq(x). 

A related process which appears for the first time at 0(a s ) is Xi production. 
Because the Yang-Landau-Pomeranchuk theorem [29| forbids gg — > Xi & t lowest (or 



any) order, the production of xi states begins at this order via the processes gg — ► 
fi'Xi) 19 ~^ QXh an d 99 ~^ QXi- While suppressed relative to X2 production by an 
additional factor of as, the relatively large radiative branching ratio to observable if) 
states (Br(xi — > ^l) = 0.27 compared to Br(x2 — ► 4>l) = 0.14) partly compensates 
so we expect Xi production to be not unreasonably small. 

Because direct if) production (via gg — > gifi) and Xi production occur at the same 
order in as, the ratio of their production cross-sections may well be less sensitive 
to assumptions about the momentum scale used in the perturbative calculations and 
so could constitute a further test of a potential model description of if) and x state 
production as well as allowing for an improved confrontation with data on VVX1/X2 
production ratios in hadroproduction experiments. The variable center-of-mass energy 
available at RHIC (y/s = 50 — 500 GeV) will be quite useful in probing the energy 
dependence of these ratios as well in probing the x-dependence of AG(x). 

In this work, we report on a set of NLO (i.e. 0(a z s )) calculations of relevance to 
these (and related) questions. We calculate a subset of the spin-independent radiative 
corrections to Xi production and include related results for xo production although 
the small radiative decay branching ratio to observable if) states make their production 
much less useful. Specifically we calculate the renormalized cross-sections for the sub- 
processes qg —* qxo,2 and qq — > gxo,2 to assess the relative importance of initial states 
containing quarks and find that their contribution is likely always smaller than 15% 



of the total cross-section. This implies that the dominant processes for \2 production 
will still be from gg initial states although a full set of NLO corrections will likely be 
necessary for a confident extraction of the polarized gluon distributions. We stress 
that we have not performed a complete set of radiative corrections to \2 production 
but we do feel that our partial results will be useful in assessing the observability of 
various x states and the likely spin-dependence of their production for the RHIC spin 
program. 

Perhaps more importantly, we also present all of the 0{a\) contributions to Xi 
production, including their spin-dependence and find that xi production increases in 
importance relative to direct if) and X2 production, providing up to 25% of the lepton 
pair yield from charmonium states at the highest RHIC energies. The spin-dependence 
of Xi production, however, is much weaker than for either direct ip and X2 production 
and may not, therefore, contribute much to measurements of AG(x). 

In the next section (2), we review the simplest ip and x production mechanisms 
and describe the calculations of the partonic level cross-sections for the Xi,2 processes 
we consider while in Sec. 3 we present numerical results for pp collisions in the energy 
range accessible to RHIC (namely 50 — 500 GeV). Finally, in Sec. 4, we discuss some 
aspects of the likely spin-dependence of these new production mechanisms and make 
some final comments. 

2 Partonic Level Cross-sections 

For reference, we begin by recalling the 0(a 2 s ) gg — > xo,2 cross-sections which form the 



basis for the leading order description of x production. They are well-known |30| and 



are given by 

<r(j99^>Xo) = Ij 7 5(1 -z) (2) 

v(gg->X2) = ^ 7 S(i-z) (3) 

where z = M 2 / s and R' P {0) is the derivative of P-state radial wavef unction at the 

origin. For completeness, we also find it useful to quote the result for gg — > gtp, 
namely 

ff(^ -» #) = 9M 5 I( < Z ) ( 4 ) 



where 



I(z) = 2z 2 



1 + z 2z\n(z) 1-- 2z 2 ln(z) 



(5) 



where -Rs(O) is the S-state radial wavefunction at the origin. (We note that, despite 
appearances, Eqn. 5 is finite as z — * 1.) 

The first of the next-to-leading order corrections we consider are those for qq — > g\j 
for J = 0,1,2. The matrix elements for these processes are readily available |K| an d 
can be integrated (over dt) to yield the finite results 

1287ra||i?^(0)| 2 2 2 (l-3^) 2 

»(*-**) - 25to ! l ^ (0)|2 ^ (1+ : ) (7) 



&(qq -> 5x2) 



27M 7 (1 - ^) 

2567ra||i^(0)| 2 ^ 2 (1 + 3z + 6^ 2 ) 



81M 7 (l-«) 

The factors of (1— z) -1 are not artifacts of any renormalization but are symptomatic 
of P-wave bound state divergences which have been noticed previously [^0| and which 



can be regulated by restricting the region of integration to s > (M + A) 2 where A is a 
typical binding energy, A ph 0.3 GeV. (Similar logarithmic divergences appear in the 



production of .D-states [3]J and in the multi-parton decays of P- and D-states [32, [33| 
as well.) We also note that the same kinematic factors appear here as in the related 
processes Z Q —>■ xj + 7 [0- New techniques for handling the related logarithms of 
binding energy which appear in P-wave decays in a rigorous way have been introduced 
J35H and could eventually prove useful in this context as well. 

The factor of (1 — 3z) 2 in Eqn. 6 implies a vanishing cross-section at some physical 
center-of-mass energy which is not obviously related to any symmetry. (Unlike radia- 
tion zeros p6 , this zero appears in the cross-section when integrated over angles. As we 



will see below, the zero may persist beyond tree level, also unlike the case of radiation 
zeros.) For future reference, we note that in each of these cases, helicity conservation 
is sufficient to determine that the partonic level asymmetries are clll — — 1- 

The next subprocesses we consider are q + g — > q + \j. The cross-sections (in four 
dimensions) for these processes have also been calculated [Q but when one attempts 



to integrate them over angles one encounters infinities, at least in the case of xo an d 
\2i which must be regulated. The Xi case is more straightforward and we treat it 



first. Because of Yang's theorem [ 29j , when the t-channel gluon (g*) becomes soft the 
amplitude describing xi ^ 99* must vanish so that there is no divergence in the phase 
space integral for this case. The (already) finite result is 

t(qg - qxi) = g^r )l [(1 -z)(5-4z- Az 2 ) - 3z 2 ln(z)}. (9) 

For the xo and xi cases we must regularize the infinities encountered in the angular 
integration and for this we use familiar dimensional regularization techniques as in Ref. 
[p^[. (A complete set of calculations for the radiative corrections to gg — > 1 S (QQ) 



have also been carried out using similar techniques f26f.) The matrix elements for qg 
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qXo,2, summed/averaged over all final/initial polarizations, must then be calculated in 
N = A — 2e dimensions. After a long calculation using FORM [B3], we find for the Xi 

case 

yi A ,|2_ 647r 2 c4|i^(0)| 2 Fi 2 \s,t,u) + eFP(s,t,u) 

Z - 1 ' ' 9M 3 (-£)(£ -M 2 ) 4 l ' 

where 

F 4 (2) (s, t, u) = {t- M 2 ) 2 {t 2 + 6M 4 ) - 2ws(t 2 - 6M 2 (£ - M 2 )) (11) 

and 

F e (2) (s,t,w) = (t - M 2 )[t 2 (M 2 - t) + 12(M 2 (M 2 - t) - us)] (12) 

while for xo we have 



327r 2 a||^(0)| 2 Ff( S) t,M)-eF(°»(s,t 



u 



| - M| " 9AP (_;)(; _ M 2 ) 4 ^ 

where 

F| 0) (s,t, M ) = (t-3M 2 )V+« 2 ) (14) 

and 

F £ (°)(s, £, u) = (t - 3M 2 ) 2 (t - M 2 ) 2 . (15) 

To be consistent with more recent treatments of radiative corrections we choose to 
average over the initial state spins by using 2(1 — e) degrees of freedom for the gluon 
which implies an additional factor of (1 — e) _1 in Eqns. 10 and 13. (This choice has 



been discussed by Ellis and Sexton |]39| and is used in Ref. pqj .) Furthermore, we 



choose to work in the MS scheme in which case the appropriate phase space factor is 

ps = h (^)^ (i - z) " 2e r dy [y{i - y)r (i6) 



where the subprocess invariants are 

M 2 ~ M 2 M 2 
s = — , t = (l-*)(l-y), u = (l-z)y. (17) 

When one performs the angular integrals, the divergent 1/e term is proportional to 
the appropriate Altarelli-Parisi splitting function (P gq (z) in this case) and is absorbed 
into the scale-dependent parton distributions leaving a finite result. Then using the 
natural choice of scale // = M we find the cross-sections 

, 167ra||i4(0)| 2 . . 
<?(qg -► QX2) = 27M^ ^ ( ] 

where 

G 2 (z) = 36(2-2z + z 2 )\n(l-z) -3z 2 ln(z) 

-53 + 69z - 18z 2 + 20z 3 (19) 



and 



where 



^(^ -> <?Xo) = 27M 7 G ° ( ^ (20) 



G (z) = 54(2 - 2z + ,2 2 ) ln(l - z) -3;z 2 ln(z) 

+4(-35 + 422-9z 2 + 2z 3 ). (21) 

We first note that, after renormalization, these contributions to the \2 an d Xo 
cross-sections are, in fact, negative which is quite similar to the case of Drell-Yan 
production |37| where the qg — > qj* correction to the tree-level process qq — > -y* is 



also negative. The contribution of the renormalized qg — » q 1 So(QQ) cross-section to 



tjq production calculated in Ref. [£6| is similarly negative. We recall that the already 
finite Xi contribution is positive. 

9 



As an aside, we note that the N = 4 — 2e matrix-element squared for qq — > g%o 
can be obtained from Eqns. 13-15 by crossing and one can see that the zero in the 
matrix element at z = 1/3 (i.e. at s — 3M 2 ) in xo production is seemingly still present 
beyond tree level. We have no explanation for this fact. 

The final case we then consider are the diagrams leading to gg — > gxj. For the 
case of xo,2 these diagrams must be combined with the virtual corrections to gg — ► ^ ,2 
to obtain a finite result and we do not attempt a complete analysis of those diagrams 
here. For the Xi case, as there is no 0{a 2 s ) contribution, the gg —> gx\ diagram gives 
the first non-vanishing contribution and it is well-behaved when integrated over angles. 
The matrix element squared for gg — * gx\ has been calculated in Ref. [4(| and the 



total cross-section can be obtained by directly integrating the four-dimensional results 
to find 

Kyy yAlJ M 7 (l + zf(l-z) A V ; 

where 

H^z) = z 9 + 39z 8 + U5z 7 + 251z 6 + ll9z 5 

-153/ - 17z 3 - U7z 2 -8^ + 10 (23) 

H[(z) = z 8 + 9z 7 + 2Qz 6 + 28z 5 + 17 z 4 

+7z 3 - 40z 2 -Az-A. (24) 

Because there is an s-channel contribution to this processes, one obtains the bound- 
state divergences (the 1/(1 — z) A factor) which are also present in the qq —>■ gxj cases 
considered earlier and which are handled in the same way. 
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3 Numerical Results 

We can now use these cross-sections to evaluate the contributions of quark (and anti- 
quark) initiated processes compared to the tree-level gluon-gluon fusion mechanism for 
Xo,2 states and all the C(a|) contributions for \i production. In all our calculations 
we use as(Q 2 = M 2 ) = 0.26 (corresponding to a leading-order A = 200 MeV) and 
the wavefunction values \R S (0) \ 2 = 0.7 GeV 3 and \R' P (0)\ 2 /M* = 0.006 GeV 3 ||. We 



use a recently updated LO parameterization of parton distributions [[|1] evaluated at 
Q 2 = M 2 and the P-wave cutoff parameter A = 0.3 GeV. The effects of changing A 

A 

by ±0.1 GeV is at most a few percent and so is not a major source of theoretical uncer- 
tainty. A dependence on A will also arise in the renormalized gg — > gxo,2 cross-sections 
and we expect a similar lack of sensitivity there as well. 

In Fig. la, we plot the contributions to \2 production in pp collisions corresponding 
to the LO gg — > \i an d NLO qg — > q\2 and qq — > g\i subprocesses. Since the 
contribution from the renormalized qg initial state is negative, its absolute value is 
plotted. One can see that the qq processes are nowhere important in the RHIC energy 
range while the qg initiated events comprise approximately 15% (30%) of the total 
(Born) level \2 cross-section (assuming a constant K = 2 factor) which is not an 
overwhelmingly large effect. The contribution from qg and qq initial states is similar 
to that found for i]q production found in Ref. |26|] although the contributions from qg 



states is somewhat larger. This implies that the total cross-section is still dominated 
by gluon fusion and so can still provide information on the polarized gluon distribution. 
A reliable extraction of AG(x), however, may well require, however, a full set of spin- 
dependent NLO corrections. 
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The various (already finite and positive) contributions to xi production from Eqns. 
7, 9, and 22, are evaluated and shown in Fig. fb. These results can then be combined 
with the direct ip production mechanism to provide estimates of the yield of e + e~ 
pairs from ip, X2, and xi states respectively once radiative branching ratios for the 
Xi,2 states are included in Fig. 2. We see that the relative amounts of if) and X2 
contributions to the observable lepton pair cross-section stays roughly constant over 
the RHIC range while the xi contributions increase to a rather substantial fraction at 
the highest RHIC energies. This fact offers yet another motivation for the usefulness of 
the high resolution electron and photon detection possible with the PHENIX detector 
at RHIC. We recall that an empirical K-factor of roughly 2 is required for the combined 
X2 and ip contributions to the observble ifj production data so that the prediction for 
Xi contribution (as well as the other two) might well be underestimated here by a 
factor of 2. 

4 Spin-spin Asymmetries 

One of the strongest motivations for considering quarkonium production at RHIC is 
the possibility of measuring the spin-dependence of the various production mecha- 
nisms, part of which includes their possible sensitivity to the longitudinally polarized 
gluon distribution. A full treatment of the spin-dependent radiative corrections to X2 
production would be necessary to fully assess the extent to which the suggestion of 
Cortes and Pire really provides a truly direct measurement of AG but, barring that, 
some useful results can be already be gleaned from our partial study. 

Since for the Xi case, all of the individual cross-sections were finite at lowest order, 
we can make direct use of the individual helicity amplitudes for the contributing pro- 
cesses to derive the partonic-level spin-spin asymmetries for all of the 0{a%) processes 
giving rise to Xi production. 
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For the case of qg — > qxi, we can integrate the individual helicity amplitudes found 



in Refs. 17 and 1H to find 



a-LL{qg -> qx\ 



3z{l - z + z\n(z)) 



(25) 



(1 - z)(5 - 4z - Az 2 ) - 3z 2 \n(z) ' 
The corresponding spin-spin asymmetry for the total gg — > g\\ cross-section is simi 
larly calculable and is given by 

z[{\ - z 2 )A x (z) + 12z\n{z)A[(z)} 



where 



a L h{gg -► gxi 



(z 2 - l)H x {z) - 12z 2 ln(z)i7i(2) 



AAz) = z 8 + 33/ + 145/ + 271z 5 + A3z 4 



+55z 3 - 273z 2 - 23z - 12 



(26) 



(27) 



A[(z) 



z 8 + 8z 7 + 25/ + 29/ + 34/ 



-32/ - 3/ - 21z - 1. 



(28) 



and the H\(z) and H[(z) are given in Eqns. 23 and 24. As mentioned above, all 
of the qq — > g\j cross-sections have an = — 1 from helicity conservation. The first 
two of these asymmetries are plotted versus the natural variable z = M 2 /s in Fig. 3 
where it can be noted that they are positive, and clearly smaller in magnitude than the 
maximally large value of cill = —1 for lowest order \2 production. For comparison, 
we also plot the corresponding asymmetry for direct ip production, namely 

" (1 - z 2 )F{z) + 2G(z) \n(z) 



a-LL 



z 2 )G(z) + 2z 2 F(z)\n(z) 



where 



F(z) = / + 2z + 5 and G{z) = / + 4/ + z + 2 



(29) 



(30) 
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which is seen to be somewhat but not dramatically larger. 

A useful measure of the average 'analyzing power' or spin-dependence in such a 

reaction is the average spin-spin asymmetry 

Stj [ dx i I dx 2 f i (x 1 , Q 2 )fj(x 2 , Q 2 )a LL da 
Yaj I dx i I dx 2 f i (x 1 , Q 2 )fj(x 2 , Q 2 )da 

which measures both the importance of the process to the total cross-section and its 
spin-spin asymmetry. (If the unpolarized parton densities, fij(x, Q 2 ), in the numerator 
were replaced by the corresponding polarized distributions, Afij(x,Q 2 ), the result 
would be the observable spin-spin asymmetry.) The average spin-spin asymmetry for 
both xi production (including all three contributing processes) as well as that for 
direct ip production are shown in Fig. 4 where we see the spin-dependence of Xi 
production is not at all large. We comment that the average value of z = M 2 /s for the 
99 ~~ * 9^ {.99 ~* 9Xi) process is roughly 0.5 (0.25 — 0.15) over the energy range plotted 
which accounts for the size and relatively constant values of < clll >. So, while the 
measurement of xi production will provide another cross check on a potential model 
description of charmonium production, it will not likely give much information on the 
polarized gluon content of the proton. 

As mentioned above, we have not evaluated the spin-dependent renormalized cross- 
sections for xi production (which in this case require a careful treatment of the helicity 
amplitudes in arbitrary dimensions, i.e. the 75 problem) but we can perhaps get 
some feel for the partonic level asymmetries for qg —>■ qxo, 2 by examining the known 
2^2 asymmetries [17| in the t — ► limit. While the partonic level cross-sections 



are divergent in this limit, necessitating the renormalization procedure above, the 
asymmetries are well-behaved. We find in this limit (when t — > and u —>■ (M 2 — s)) 

—2z + z 2 

d L L{qg -» QX2) = - — - — ; — = -dLiiqg -> qXo)- (32) 

2 — 2z + z z 

14 



where z = M 2 /s. We plot this function in Fig. 5 and note that at threshold (z = 1) 
the partonic level asymmetry is equal to — 1(+1) just as for gg — > X2(Xo) & t tree level. 
For comparison, we also plot on the same figure the partonic level asymmetry for 
99 ~~* 9X2, again in the same limit, and note a similar behavior. As before, the Xo,2 
asymmetries are trivially related and we find 

, , v -z(2 -3z + 2z 2 ) A 

aLL{99 -» 9X2) = —7i , 2 , 2 = -aLL{99 -► 9Xo)- (33) 

(i — z + z i y 

The purely gluonic asymmetry, however, stays near its threshold value of —1 over 
a wider range of z and so will contribute more to the observed asymmetry. (These 
expressions are also found in the u — > limit where divergences also occur in the 
99 ~> 9Xo,2 processes.) If the fully renormalized gluon induced cross-section were 
described by these limiting values for the asymmetries, then an average value of < z > 
of 0.5 (0.2), as for gg —> gip {gg — > gxi), would correspond to < clll >~ —0.9 (—0.5). 
These results suggest that the spin-dependence of the NLO gg and qg processes will 
be somewhat 'softer' but will still retain much of the same general structure (most 
especially the sign) of the tree-level spin-spin asymmetry in contrast to, for example, 



open heavy flavor production as discussed in Ref. |L2 . 

As a final comment, we note that X2{cc) production RHIC probes the polarized 
gluon distributions in a kinematic region given by 

VT= —= ^0.07- 0.007 . 

We note that the production of a 1 S (tt) toponium state (presumably detected via 
its two-photon decay) also has a maximally large partonic level spin-spin asymmetry 
given by clll = +1- A future supercollider (such as the LHC) with a polarization 
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option would be expected to probe r\ t production in a kinematic range given by \fr ~ 
.0017 — .0025 which would then be in the region of polarized gluon densities already 
probed by RHIC charmonium studies. Since the value of as probed by such heavy mass 
states would be much smaller than for charmonium, the effects of radiative corrections 
would be expected to be much smaller and the interplay between polarized gluon 
densities and the leading-order production mechanisms would be expected to be much 
more direct. In the same context, standard model Higgs boson production via gluon 
fusion, i.e. gg — > H via quark loops, has the same leading-order spin-spin asymmetry, 
namely +1. Higgs boson production at the LHC via this mechanism would then be 
in the same kinematic region as X2(cc) production at RHIC for Higgs boson masses in 
the range 200 — 2000 GeV which covers most of the expected range. In this regard, the 



spin-independent radiative correction calculations of Refs. j|3], |44| are already useful 
in assessing the contributions of various NLO processes. 
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Figure Captions 

Fig. 1 The total cross-section, a(nb), for X2,i (cc) production in pp collisions versus y/s 
(GeV). For 1(a), we plot the results for X2 production with the solid (dotted, 
dashed) curves corresponding to the 0(a 2 s ) gg — ► X2 (qg —> qX2, qq -^ 9X2) 
contributions. The negative of the qg contribution is plotted for simplicity. For 
1(b), the results for xi production are shown where the solid (dotted, dashed) 
curves correspond to the gg — > gx\ (qg —> qXi-, qq ~ *■ 9Xi) contributions. 

Fig. 2. The total cross-section for i[) production (times leptonic branching ratio, Br(ip — ► 
e + e~)) in pp collisions versus y/s(GeV). The contributions from the lowest order 
99 -^ 9^ (solid), 0(a 2 s ) gg — ► X2 (dotted), and total 0(a%) Xi contributions 
(dashed) from g, qg, and qq initial states are plotted. 

Fig. 3. The partonic level spin-spin asymmetry, a LL , (in the total integrated cross- 
section) versus z = M 2 / 'y/s for xi production. The solid (dotted) curves corre- 
spond to the asymmetry for gg — ► gx\ {qg —> qXi)- The asymmetry for direct 
ip production via gg — > gip production (dashed curve) is also shown for compar- 
ison. The partonic level asymmetries for qq — > (7X0,1,2 and for gg — > X2 are both 
aLL = -1- 

Fig. 4. The average spin-spin asymmetry, < (lll > versus the center of mass energy, y/s 
GeV for pp collisions. The results for total Xi production (dots) (from gg, qg, 
and qq processes) and direct ip production (solid) (via gg — ► gip) are shown. 

Fig. 5. The partonic level spin-spin asymmetries, o,ll, versus z = M 2 /s for X2 produc- 
tion in the limit that t — ► 0. The solid (dotted) curves correspond to the processes 

99 -»■ 9x2 (qg -> 9x2)- 
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